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Abstract 

We study the transverse momentum (pt) spectrum of charged particles produced 
in deep inelastic scattering (DIS) at small Bjorken x in the central region between the 
current jet and the proton remnants. We calculate the spectrum at large pT with the 
BFKL ln(l/x) resummation included and then repeat the calculation with it omitted. 
We find that data favour the former. We normalize our BFKL predictions by comparing 
with HERA data for DIS containing a forward jet. The shape of the x distribution of DIS 
+ jet data are also well described by BFKL dynamics. 



^On leave from H. Niewodniczanski Institute of Nuclear Physics, Department of Theoretical Physics, ul. 
Radzikowskiego 152, 31-342 Krakow, Poland. 



1. Introduction 



An intriguing feature of the measurements at HERA in the small x domain is the possible 
existence of significant ln(l/x) effects. A major part of the rise observed for the structure 
function F2 with decreasing x may be attributed to the resummation of the leading ln(l/x) 
'BFKL' contributions. An excellent unified BFKL/GLAP fit of F2 in the HERA regime has 
recently been obtained using a "flat in x" input and the rise due to BFKL- type effects has 
been quantified within this description. However, the growth of F2 with decreasing x can be 
described equally well by pure GLAP ||^ In(Q^) evolution from suitably chosen input parton 
distributions so the main origin of the rise is still an open question. The observable F2 is too 
inclusive to distinguish between these alternatives. The study of deep inelastic scattering (DIS) 
events containing an isolated forward jet [|, |^ is a better discriminator of the underlying small 
X dynamics. The process is sketched in Fig. 1(a). In this case we effectively study DIS off known 
parton distributions and so we avoid the ambiguity in the choice of the input distributions. The 
method is theoretically attractive. The summation of the leading ln(l/x) contributions gives an 
{x/xj)~^ behaviour of the BFKL ladder connecting the photon to parton a. Here x is Bjorken 
X and Xj is the fraction of the proton's longitudinal momentum carried by the parton jet. An 
unambiguous measurement of the exponent A looks feasible. In practice a major problem is 
the identification of the jet due to parton a, and the measurement of its momentum, when it is 
close to the remnants of the proton. Typically the clean observation of the jet requires Xj ^0.1 
and so in this process we lose about a factor of 10 in the 'small x reach' of HERA. 

Besides the x~'^ growth as x decreases along the BFKL ladder, a second characteristic 
feature is the diffusion in In where kx are the transverse momenta of the gluons emitted 
along the chain. One way the diffusion manifests itself is in an enhancement of the transverse 
energy (Ft) flow in the central region between the current jet and the proton remnants see 
Fig. 1(b). In principle the diffusion can enhance Ft from both the 'upper' and 'lower' BFKL 
gluon ladders, which are denoted by $ and / in Fig. 1(b). However, the x reach at HERA is 
insufficient to fully develop the In k^ diffusion in both ladders simultaneously. Nevertheless, the 
effect is quite appreciable giving at the parton level an energy flow F^ < 2 GeV/unit of rapidity. 
However the clean parton level prediction can in practice be masked or mimicked by the effects 
of hadronization. Thus, although the prediction for Ft is in agreement with observations 0] 
we cannot definitely conclude that it is due to ln(l/x) resummations. 

An interesting way to overcome this ambiguity is to consider the emission of single parti- 
cles at relatively large transverse momentum pt in the central region The single particle 
spectrum at sufficiently large values of pt should be much more immune from hadronization 
and more directly reflect the In k^ diffusion from the BFKL ladders. 

The outline of the contents of the paper is as follows. In Sec. 2 we use the data for the 
process DIS + forward jet to normalise the BFKL function $ shown in Fig. 1(a). To be precise 
we numerically solve the BFKL equation for $ using the amplitude for the quark box 
(and crossed box) as input at a value zq oi z = x/xj which is chosen so that the resulting $ 
reproduces the DIS -|- jet data. Also, for completeness, we present in Sec. 2, an analytic form 
for $ which is valid for fixed as-, and which has been the basis for a recent analysis. In Sec. 3 
we give the formula necessary to calculate the transverse momentum (pj-) spectrum of single 
particles produced in the central region. The process is shown in Fig. 1(c). The predictions for 
the Pt spectra (with and without the BFKL effects included) are compared with HERA data. 
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Finally in Sec. 4 we give our conclusions. 



2. DIS + forward jet events 

We first calculate the cross section for DIS containing a forward identified jet. This so-called 
"Mueller" process is a valuable probe of small x dynamics in its own right. We compare with 
HERA data to normalise the function $ shown in Fig. 1(a). There are uncertainties in the 
normalisation, and even the shape of the x distribution is dependent on subleading ln(l/x) 
corrections. 

The variables of the process are shown in Fig. 2. As usual the variables x and y are given 
hj X = Q'^/2p ■ q and y = p ■ q/pe ■ P where p, Pe and q denote the four momenta of the proton, 
the incident electron and the virtual photon respectively, and = —q^. The variables Xj and 
kjT are the longitudinal momentum fraction and transverse momentum carried by the forward 
jet. The differential cross section is given by |Q 



- Ux, fde 



dxdQ^ J ^'^^ J xQ^ 



dF2 1 2 dFr 



dxjdkjj, 2 dxjdk'jrp 



(1) 



where the differential structure functions have the following form 

mr-^^f^(^,^^^-{h^r:<A (2) 



for i = T,L and F2 = Ft + F^. We have assumed strong ordering at the parton a - gluon 
vertex. Assuming also t-channel pole dominance the sum over the parton distributions is given 
by 

E/a = ^? + |E(? + ?1- (3) 

a <] 

Recall that these parton distributions are to be evaluated at {xj, kjj.) where they are well-known 
from the global analyses, so there are no ambiguities arising from a non-perturbative input. 



The functions ^i{x/xj, k'jrp, Q^) describe the virtual 7 + virtual gluon fusion process includ- 
ing the ladder formed from the gluon chain of Fig. 2. They can be obtained by solving the 
BFKL equations 

$,(z,4,Q^) = <l>fH^,4,Q2)+ 

«5 /' ^ / ^ [HA {q + M^ Q')) - HA 4, Q')0(4 - q')] (4) 

J z Z J TTq 

where as = 3^5/71. The inhomogeneous or driving terms correspond to the sum of the 
quark box and crossed-box contributions. For small z we have 

^f\z,klQ') ^ $r(z = o,4,Q^) - $r(4,Q^). (5) 

We evaluate the ^f'^ by expanding the four momentum in terms of the basic light-like four 
momenta p and q' = q + xp. For example, the quark momentum k, in the box (see Fig. 2) has 
the Sudakov decomposition 

K = ap — Pq' + kt- 
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We carry out the integration over the box diagrams, subject to the quark mass-shell constraints, 
and find 



2E 



^2 «g Q'^ 
+ (1 - (if 



dp J d^KT 



Pi 



D^2 



2E 



.2 

q " 'T 



^« 7r2 A;2 7o 
where the denominators Dj are of the form 



dp / ^2^^/32(1-/3)^ 



= 4 + /3(i - /?) q2 + 



+ m 
' 1 



1 



1 



1 \ 



(6) 



(7) 



The light m, and s quarks are taken to be massless (m^ = 0) and the charm quark to have 
mass irir = 1.4 GeV. 



2.1 Analytic form of $ for fixed 0:5^ 

We solve the BFKL equation for $ numerically, which allows the use of running as and the 
inclusion of a charm quark mass. However, it is informative to recall the analytic solution which 
can be obtained if as is fixed and we assume that the quarks are massless. The first step is to 
rewrite the driving terms for = in the form 



\ , j/^^ + (1 - /^n [A^ + (i-A)2] 

^ [A(1-A)A;2 +/3(1-/?)Q2] 



A(l-A)/3(l-/3) 



where A is the Feynman parameter which appears in the representation 



o''^[ad7 

.(0) 



(8) 
(9) 

(10) 



We see that, for fixed 0:5 and niq = 0, the $• are functions of a single dimensionless variable 
r = Q'^/k'^. We may therefore represent the driving terms ^f'\Q'^/k'^) in terms of their Mellin 



transforms ^f'\'y) 



$(o)fr) 



27ri Ji-ioo 



d7 ^'\l)r 



(11) 



where i = L,T and r = Q'^/k^. The Mellin transform is useful since it diagonalizes the BFKL 
equation (^). The solutions for fixed coupling as may therefore be written 



$,(z,4,q2) 



27ri Ji-ioo 



2\ 7 



rf7(^) exp{asK{j)\n-J^f'\j) 



(12) 



where as = Sas/ir and -^'(7) is the Melhii transform of the kernel of the BFKL equation 

K{-^) = 2^(1) -^(7) -^(1-7) (13) 

with \I/(7) = r'(7)/r(7). The functions ^f'\'y) are obtained by inserting (|]) and @ into the 
inverse relation to (|TT|). We find 



4°^ (7) 



as 



^ 47r Jo Jo Jo [A(l - A) + - p)r] 

2 (^s 
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(14) 
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''^^"Vo'^Vo'^^[A(l-A)+/5(l-/3)r] 



5(-7 + 2,-7 + 2)5(7 + l, 7 + 1: 



(15) 



where B{x,y) = r{x)r{y)/r{x + y). The derivation of the analytic formula relies on as being 
fixed. This approach has been used by Bartels et al. 0] to estimate the DIS + forward jet cross 
section taking the coupling as{k^) in formulae ([T2|). The prediction has the general shape of 
the HI data as a function of x, but the calculated cross section exceeds the data by some 20% 



In the z limit the formulae reduce to the conventional z ^ BFKL behaviour 



$f(7=|) 



k^J {asK"a)\nl/zy2 



(16) 



where for simplicity we have omitted the Gaussian diffusion factor in ln(A;|./(5^). If we evaluate 
the various functions at 7 = | we obtain 



.-A 



512 ^21C(3)/2U^/ ^Hl/z) 



1 + 



.In(lA), 
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(17) 



where \ = as K{^) = 4 In 2. 



2.2 Normalisation of $ and the description of DIS + jet data 

Our calculation of the DIS + forward jet process differs from that of ref. |p in that we numeri- 
cally solve the BFKL equations. Therefore, we are able to explicitly include the rric 7^ charm 
contribution. We also allow the coupling 0:5 to run. To be precise we solve the BFKL equation 
(§) rewritten in terms of the modified function as{k'^T)^ii^y ^t; Q^) following the prescription 
that was used in ref. This choice of scale for as is consistent with the double logarithm 



limit and with the NLO ln(l/x) analysis of ref. 
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We determine the functions $j for z < Zohy solving the BFKL equation as described in 
starting from the boundary condition 

n^o,k'r,Q') = $r(^o,4,Q') ^ (18) 

where ^f^\k^,Q'^) are the contributions of the quark box (and crossed box) given in (^. We 
take u, d, s to be massless and the charm quark to have mass rric = 1.4GeV in the summation 
over the quarks. We then use ([|) integrated over x and and (0) to calculate the DIS + 
forward jet rate corresponding to the cuts used in the HI measurement. That is the forward 
jet is constrained to the region 

7° < Oj < 20°, Ej > 28.7GeV, kjT > 3.5GeV, 

whereas the outgoing electron must lie in the domain 

160° <e'^< 173°, E'^ > llGeV, y > 0.1 

in the HERA frame. Finally HI require ^Q^ < A;|j. < 2Q^. The BFKL calculation is compared 
with the data ||10| in bins of size Ax = 5 x 10~^ in Fig. 3. The parameter zq is adjusted to give 



a satisfactory normalization of the calculation. We find that the HI data require zq = 0.15. 
The predicted shape of the distribution is in good agreement with the data. 



3. Single particle pt spectra 

We first use Fig. 1(c) to obtain the differential cross section for the production of a hadron 
of transverse momentum px and longitudinal momentum fraction Xh- Then we calculate the 



charged particle spectra relevant to the recent observations at HERA [12 



3.1 The cross section for charged particle production 

The cross-section for single particle production is obtained by convoluting the inclusive 
cross-section for the production of a single parton with the parton fragmentation function. 
The differential cross section for the inclusive production of a single parton of longitudinal 
momentum fraction xj and transverse momentum kjT has the generic form of (|l]). We have 

dxj dkjrp dx dQ'^ xQ'^ 

Now for small x, and in the central region away from the current jet and the proton remnants, 
we expect gluonic partons to dominate where the gluons are radiated within the BFKL ladder. 
The differential structure functions occurring in (|19]) are then given by 

dFi r d?hp r d?k, 

dx, dk] ^ J ^ J 1^ 

with i = T,L and where for simplicity we have omitted the subscript T from the gluon transverse 
momenta, kjT, kpT and k^T, see Fig. 4. The functions $j are those of Sec. 2 which control the 



. dF2 ^ ly2_9FT_ 
dxj dk'^rp 2 dxj dk^j^ 



(19) 



as{k^)k'pk 



fi.Xj,kl) $i 



X 



ktQ^\5\k, 



"7' 



kip kf-y ^ 



(20) 
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DIS + forward jet rate, while / is the unintegrated gluon distribution which satisfies the BFKL 
equation 



— z 



dz 



as 



(21) 



The expression in square brackets in (|20|) arises from the (square of the) BFKL vertex for real 
gluon emission, see Fig. 4. 

In practice we evolve ( pi]) down in z from the boundary condition 



fizo,k') = r^(zo,^') 



d 



zo g^^{zo, A;2) 



d\n{kyki) 



(22) 



where here z = zq with zq = 10 ^, and where g^^ is the conventional gluon distribution 
obtained from a global set of partons. As before we allow the coupling to run, that is we take 



m 



Moreover, we impose an infrared cut-off /cq = IGeV . That is we require the 



arguments of / to satisfy k > k^ and {k + q) > k^. Similarly, the integrations in ( PO] ) are 
restricted to the regions fc^, k"^ > k^. We may include the contribution AFi from the region 
kp < /cq by assuming the strong ordering approximation, k'i <^ k'i ^ k'^- , at the gluon vertex. 



This contribution to (pOD then becomes 

d{AFi) 



X 



^ dxjdkj 



1^2 Jl.2 / 



kj 



X 



kj , Q 



(23) 



Most of the time, however, for the calculation relevant to the HERA data, the variable Xj is 
not small enough for the BFKL equation to be applicable for the function /. In these cases, 
that is when Xj > zq, we therefore again assume strong ordering k"^ <^ k"^ r-^ k"^- . In addition we 
include the contributions from quark and antiquark jets. We then obtain 



dxjdkj 



asm ^ 



k] 



X 



, kj, Q 



(24) 



where the parton distributions are to be evaluated at {xjjk"^) The differential cross section 
for single particle (h) production is obtained by convoluting the jet cross section with the 
fragmentation functions D for the parton h transition 



dxhdp^dxdQ'^ 



L I "^""^ I ^^""^ " '"^'^^^^^ " '^'^ { dx/k]dxdQ^ ""'^^'^'^ 



+ 



9 ^ 

» q 



dxjdkWxdQ"^ 



D'{z,f.') + 



dag 



dxjdkjdxdQ' 



(25) 



where ag and cxg are the contributions to the cross section aj for gluon and quark and antiquark 
jets respectively. The fragmentation scale //^ is of the order of /c? The cross section for charged 
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particle production is obtained by summing over all possible charged hadrons h. 



3.2 Predictions for the single particle pt spectra 

The data for the single (charged) particle pt spectra are presented in the form {dn/dpT)/N 
where n is the multiplicity and the total number of charged particles in a given x, bin 



10|. To calculate this pt spectrum we evaluate 



N dpT ~ \^ dpT dx (9g2 ) I dxdQ^' ^ ^ 

where dah/dpTdxdQ"^ is obtained from ( ^5] ) by integrating over Xh- We take the central values 
of x, in the bin. The integration limits are fixed by the limits on the pseudorapidity interval 
under consideration. To be precise we use 



= Jt^Ptc-^ (27) 




where t] is the pseudorapidity of the charged particle, i] = — lntan(6'/2) with 6 the angle with 
respect to the virtual photon direction. Finally we calculate the total differential cross section 
datot / dxdQ'^ in ( p6D from the structure functions F2 and Fl given by the MRS(R2) ||13[ set of 
parton distributions. 

Our aim is to make an absolute BFKL-based prediction to compare with the pt spectra ob- 
served by the HI collaboration. There is, however, an inherent uncertainty in the normalisation 
due to the imposition of an infrared cut-off on the BFKL transverse momentum integrations (or 
due to other possible treatments of the non-perturbative region). To overcome this problem we 
follow the procedure described in Sec. 2.2 and fix the parameters occurring in the calculation 
of the BFKL functions by requiring the prediction for the DIS + jet cross section to give the 
correct normalization of the HI forward jet measurements. The next step is to use the functions 



$j obtained in this way in the computation of the differential structure functions from (pOD 
(p3|) and (0). In this way we are able to calculate a normalized pt spectrum from ([26|). 



The BFKL prediction for the single particle spectra may be compared with the result which 
would be obtained if the BFKL gluon radiation is neglected. That is in (pO|), ( p3|) and (12^) we 
replace the functions $j which describe the solution of the BFKL equation with the boundary 
condition given by the quark box with the quark box ^f^^ only. In addition we now also assume 
strong ordering for Xj < Zo and carry out the kp integration in (pO]). This amounts to assuming 
that in a fixed-order treatment the dominant subprocess is 7(7 qqg. In our calculation the 
K integration is infrared finite since we allow for the virtuality of the incoming and exchanged 
gluons. 

So we are now in the position to give a BFKL prediction for the single particle spectra 
which can be compared with the HI data. In their measurement the HI collaboration collected 
data in nine different kinematic bins in two pseudorapidity intervals. We will focus on the three 
smallest x bins where BFKL effects should become visible. Also we will only show results for 
the lower pseudorapidity interval, 0.5 < r] < 1.5, where we expect no contamination due to 
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the fragmentation of the current jet which has not been included in the calculation. In the 
computation of the px spectra we use 



= 27.5GeV, Ep = 820GeV 

and impose the cuts which where used in the HI measurement, i.e. we require the outgoing 
electron to lie in the region 

157° <9'^< 173°, E'^ > 12GeV, y > 0.05 

in the HERA frame. Also we subtract 10% off the total cross section atot to account for 
diffractive events with large rapidity gaps which have been excluded from the measurement. 
Finally, in the sum over the charged hadrons h in (^6]) we include vr^ and K^, and we use 



the next-to- leading order fragmentation functions by Binnewies et al. |]T^. In Fig. 5 we show 
predictions for the charged particle pt spectrum in kinematic bin 1 of the HI analysis with 
central values x = 1.6 x 10~^ and = 7 GeV^. We compare the results when BFKL small 
X resummation is included in the calculation with the case when gluon radiation is neglected. 
In both cases we demonstrate the effect of changing the fragmentation scale from fi"^ = A;| to 
/x^ = {2k jY. We see that the BFKL prediction gives a good description of both the shape and 
the normalization of the HI data. On the other hand, when the BFKL effects are neglected 
the predictions lie considerably below the data. Also we see, as expected, that the spectrum 
decreases more rapidly with pt than when the BFKL resummation is included. For example 
for Pt = 1.5 GeV the two predictions differ by a factor 3.6, whereas for pt = Q GeV this factor 
is almost 10. This is a reflection of the diffusion in Ink^ along the BFKL ladder. 

The same general behaviour is seen in Figs. 6 and 7 where we show the comparison for 
kinematic bins 2 and 3, with central values x = 2.9 x 10~^, = 9 GeV^ and x = 3.7 x 10~^, 

= 13 GeV^, respectively. We find that in all three small x bins of the HI analysis the 
data strongly support the inclusion of BFKL resummation in the calculation of the pt spectra. 
Reasonable variations of the fragmentation scale do not allow for a description of the data when 
BFKL effects are neglected. 



4. Conclusion 



We studied the DIS + forward jet process including massive charm in the quark box and 
solving the BFKL equation numerically for running coupling. We found that BFKL dynamics 
describe the shape of the x distribution of the HERA data well. Next we used these data to fix 
the normalization of the solution of the BFKL equation with the boundary condition given by 
the quark box. This enabled us to give an absolute prediction for charged particle transverse 
momentum spectra at small x. We calculated the spectrum for large values of pt first including 
BFKL small x resummation in the calculation and second neglecting gluon radiation. It turned 
out that the BFKL prediction agrees well with the HI data both in shape and normalization, 
whereas the approximate fixed order result underestimates the data and decreases too rapidly 

^Even though we have normahsed $ to the DIS + jet data, there still remains some residual uncertainty in 
the overall normalisation associated with the choice of infrared cut-off used in the kj integration in (po|). Our 
results are shown for the natural choice fc^n = 1 GeV^. 
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with pt- We therefore conclude that we found evidence for the existence of ln(l/x) effects and 
for the diffusion in In which accompanies BFKL evolution. Despite these encouraging results 
it would, however, still be useful to compare the BFKL prediction for the spectrum with 
the result of the complete fixed order calculation. Experimental data for higher values of pt 
would allow an even clearer distinction between the different predictions. BFKL effects would 
also become more apparent in the pseudorapidity interval —0.5 < rj < 0.5 which corresponds 
to higher values of Xj and therefore to a longer BFKL evolution starting from the quark box. 
Of course higher Xj also means less BFKL evolution from the proton end. This is, however, 
not a disadvantage, since already for the pseudorapidity interval which we considered the main 
contribution to the spectrum comes from the region Xj > zq. We conclude that although more 
experimental data especially for higher values of px would be useful, the existing spectra show 
the presence of BFKL effects at small x at HERA. 
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Figure Captions 



Fig. 1 Diagrammatic representation of (a) the deep inelastic + forward jet, (b) the Ex flow, and 
(c) the single particle spectrum measurement. 

Fig. 2 Diagrammatic representation of a deep inelastic + forward jet event. 



Fig. 3 The deep inelastic + forward jet cross section in pb integrated over bins of size 5 x 10 



-4 



in X compared to the HI data presented at the Warsaw conference |10|. As in the HI 
measurement the forward jet was required to fulfil 7° < 9j < 20°, Ej > 28.7 GeV, and 
kjT > 3.5 GeV. The electron acceptance region is limited by 160° < 0'^ < 173°, E'^ > 11 
GeV, and y > 0.1 in the HERA frame. 

Fig. 4 Diagrammatic representation of the cross section for emission of a high transverse mo- 
mentum prp particle. 

Fig. 5 The transverse momentum spectrum of charged particles {7c~^ ,tt~ , , K~) in the pseu- 
dorapidity interval 0.5<r7<1.5in the virtual photon-proton centre-of-mass frame. The 
results are shown for kinematic bin 1 with the central values x = 1.6 x 10~^ and = 7 
GeV^. The continuous and the dashed curve show the spectra obtained with $j and / 
calculated from the BFKL equation. They only differ in the choice of fragmentation scale: 
for the continuous curve the fragmentation functions were evaluated at scale /i^ = (2/cj)^ 
and for the dashed curve at scale fi^ = kj. When BFKL radiation is neglected in the 

calculation of the pt spectra, i.e. when the quark box approximation $j = ^f'^ is used 
and strong ordering at the gluon vertex is assumed, then the dash-dotted and dotted 
curves are obtained. The fragmentation functions were evaluated at scales 2kj and kj re- 
spectively. The data points shown are from the HI measurement of the charged particle 
spectra |1T2| . 



Fig. 6 As Fig. 5, but for kinematic bin 2, x = 2.9 x 10 and Q'^ = 9 GeV^. 
Fig. 7 As Fig. 5, but for kinematic bin 3, x = 3.7 x 10"^ and = 13 GeV^. 
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